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ABSTRACT 

We here investigate the possibihty that the ultra-high energy cosmic ray 
(UHECR) events observed above the GZK hmit are mostly protons accelerated 
in reconnection sites just above the magnetosphere of newborn millisecond 
pulsars which are originated by accretion induced collapse (AlC). We formulate 
the requirements for the acceleration mechanism and show that AlC-pulsars 
with surface magnetic fields 10^^ G < -B* ^ 10^^ G and spin periods 1 ms 
^ < 60 ms, are able to accelerate particles to energies > 10^° eV. Because 
the expected rate of AIC sources in our Galaxy is very small (~ 10^^ yr^^), the 
corresponding contribution to the flux of UHECRs is neglegible, and the total 
flux is given by the integrated contribution from AIC sources produced by the 
distribution of galaxies located within the distance which is unaffected by the 
GZK cutoff (~ 50 Mpc). We find that reconnection should convert a fraction 
^ ^ 0.1 of magnetic energy into UHECR in order to reproduce the observed 
flux. 

Subject headings: acceleration of particles - stars: magnetic fields - 
pulsars: general- stars: neutron - white dwarfs 
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1. Introduction 

The detection of cosmic ray events with energies beyond 10^°eV by AGASA (Takeda 
et al. 1999), Fly's Eye (Bird et al. 1995), and Haverah Park (Lawrence, Reid, & Watson 
1991) experiments still poses a challenge for the understanding of their nature and sources. 
These ultra-high energy cosmic rays (UHECRs) show no major differences in their air 
shower characteristics to cosmic rays at lower energies and thus one would expect them to 
be mostly protons (Protheroe 1999). If UHECRs arc charged particles, or protons, then 
they should be affected by the expected Greisen-Zatsepin-Kuzmin (GZK) energy cutoff 
(~ 5 X 10^^ eV), which is due to photomeson production by interactions with the cosmic 
microwave background radiation, unless they are originated at distances closer than about 
50 Mpc (e.g., Protheroe & Johnson 1995, Medina Tanco, de Gouveia Dal Pino & Horvath 
1997). On the other hand, if the UHECRs are mostly protons from nearby sources (located 
within ~ 50 Mpc), then the arrival directions of the events should point toward their 
sources since they arc expected to be little deflected by the intergalactic and Galactic 
magnetic fields (e.g., Stanev 1997, Medina Tanco, de Gouveia Dal Pino & Horvath 1998). 
The present data shows no significant large-scale anisotropy in the distribution related to 
the Galactic disk or the local distribution of galaxies, although some clusters of events seem 
to point to the supergalactic plane (Takeda et al. 1999, Medina Tanco 1998). 

A number of source candidates and acceleration mechanisms have been invoked but 
all of them have their shortcomings (see, e.g., Olinto 2000 for a review). Particles can, 
in principle, extract the required energies from an induced e.m.f. in a circuit connected 
between the polar and the last open field line of a rapidly rotating pulsar, although it is 
not clear how the large voltages can be maintained (e.g., Hillas 1984), or be accelerated in 
reconnection sites of magnetic loops, if these can be produced, e.g., by Parker instability, 
on the surface of a pulsar (Medina Tanco, de Gouveia Dal Pino & Horvath 1997), but the 
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accelerated particles will probably lose most of their energy gain by curvature radiation 
while dragged along by the magnetic dipole field (Sorrell 1987). Alternatively, Olinto, 
Epstein, & Blasi (1999) have recently proposed that UHECRs could be iron nuclei stripped 
by strong electric fields from the surface of highly magnetized neutron stars and accelerated 
in a relativistic MHD wind. It is not clear however, how the accelerated particles can escape 
from the magnetosphere of the star, for although efficient power extraction may be possible, 
there is a dense positron-electron plasma generated with the relativistic wind (Gallant & 
Arons 1994) that will possibly modify the electric fields and also interact with the energetic 
particles. Their model also predicts that a correlation with the Galactic plane should 
become evident as data collection at the highest energies improves. In this paper we discuss 
an alternative model in which UHECRs are accelerated in magnetic reconnection sites 
outside the magnetosphere of very young millisecond pulsars being produced by accretion 
induced collapse (AIC) of a white dwarf. 

2. The model 

When a white dwarf reaches the critical Chandrasekhar mass ~ 1.4 Mq through 
mass accretion, in some cases it does not explode into a type la supernova, but instead 
collapses directly to a neutron star (e.g., Woosley & Baron 1992 and references therein). 
The accretion flow spins up the star and conflncs the magnetosphere to a radius Rx where 
plasma stress in the accretion disk and magnetic stress balance (Arons 1993). At this 
radius, which also defines the inner radius of the accretion disk, the equatorial fiow will 
divert into a funnel infiow along the closed field-lines toward the star (Gosh and Lamb 
1978), and a centrifugally driven wind outflow (Arons 1986). Recently, Shu et al. (1994, 
1999) have studied the detailed field geometry of magnetized stars accreting matter from 
a disk. Two surfaces of null poloidal field lines are required to mediate the geometry of 
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dipole-like field lines of the star with those opened by the wind and those trapped by 
the funnel inflow emanating from the Rx region. Labeled as ^''helmet streamer^'' and 
" reconnection ring" in Figure 1, these magnetic null surfaces begin or end on Y points. 
Across each null surface, the poloidal field suffers a sharp reversal of direction. Dissipation 
of the large electric currents that develop along the null surfaces will lead to reconnection 
of the oppositely directed field lines (e.g., Biskamp 1997, Lazarian & Vishniac 1999). 
Helmet streamers (or flare loops) are also present in the magnetic field configuration of 
the solar corona. The magnetic energy released by reconnection in the helmet streamer 
drives violent outward motions in the surrounding plasma that accelerate copious amounts 
of solar cosmic rays without producing many photons (Reames 1995). A similar process 
may take place in the helmet streamer of young born AlC-pulsars and the magnetic energy 
released may accelerate particles to the UHEs. The reconnection rings in the disk would 
also, in principle, be able to accelerate UHECRs, however, as we will see below, the intense 
radiation field produced in the disk must prevent accelerated UHE particles to escape from 
it. The particular mechanism of particle acceleration during the reconnection events is 
still unclear in spite of numerous attempts to solve the problem (see LaRosa et al 1996, 
Litvinenko 1996). Cosmic rays from the Sun confirm that the process is sufficiently efficient 
in spite of the apparent theoretical difficulties for its explanation. Solar flares observations 
indicate that the reconnection speed can be as high as one tenth of the Alfven velocity. As 
in the conditions we deal with this speed approaches c, the expected acceleration rate is 
large. We discuss the details of particle acceleration during reconnection events in Lazarian 
& de Gouveia Dal Pino (2000). We argue that protons can be accelerated by the large 
induced electric field within the reconnection region (e.g., Haswell, Tajima, & Sakai 1992, 
Litvinenko 1996) over a time scale ~ ARx/c, where ARx is the size of the reconnection 
zone (see below). The required electric field is less than the critical value for pair production 
and therefore is sustainable (see also §3). 
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For a Keplerian disk, the inner disk edge Rx rotates at an angular speed {GM^/R^Y^'^, 



and equilibrium between gravity and centrifugal force at Rx will lead to co-rotation of the 
star with the inner disk edge, i.e., Rx = {GMi,/Qiy^^ , which for typical millisecond pulsars 
with rotation periods = 27r/Qyt ~ 1.5 — 10 ms, mass M^, ~ 1 M©, and radius i?* = 10^ 
cm, gives Rx — (2 to 7) x 10^ cm Rq, where Rq = -R*/10^ cm. 

The primary condition usually assumed on the region to accelerate particles of 
charge Ze to energies E is that its width ARx > 2rL, where is the particle Larmour 

radius = E/ZeBx (Hillas 1984) and Bx is the magnetic field (normal to particle 

1/2 

velocity) at the Rx region, Bx ^ Bdipoie{Rx) (^7) (e-g-, Arons 1993), where 
Bdipoie{Rx) — -B* {R^/RxY is the magnetic field that would be present in the absence of the 
shielding disk, and is the magnetic field at the surface of the star. While this condition 
on Ai?x is usually invoked to allow particles to bounce back and forth thus gaining energy, 
we find that for accelerated protons the synchrotron losses may be too large if they bounce 
within the reconnection zone. Then, in our model the condition above is invoked to assure 
that the field Bx will focus particles to move within a small angle into the reconnection 
zone. Besides, one should also expect that: ARx/Rx « 1- Both conditions above imply 



a given ratio ARx / Rx and particle energy E, the stellar magnetic field B^ must satisfy 



where we have assumed = IMq, R^ = Rq, E20 = ^/lO^^ eV, ^2.5^ = 0^/2.5 x 10^ s'^, 
and B12, = 5^/10^^ G. The corresponding allowed zones in the -B* — plane are shown 
in Figure 2 for £^20 = 1 and £^20 = 10, and different values of the ratio ARx/Rx- The 
curves with ARx / Rx — 1 determine extreme lower bounds on the stellar surface magnetic 
field and the angular speed. We note that stellar magnetic fields 10^^ G < B^ <^ 10^^ G 
and angular speeds 4 x 10^ s~^ ^fl^> 10^ s~^, are able to accelerate particles to energies 



that 1 >> 



(^) > 4e-2 (GM^y E^ fi-^/l This relation indicates that for 




(1) 
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^20 ^ 1- The values above are perfectly compatible with the parameters of young pulsars 
and Eq. (1) is thus a good representation of the typical conditions required for particle 
acceleration to the UHEs in reconnection zones of AIC- pulsars. The substitution of Eq. (1) 
into the equation for Bx imphes a magnetic field in the acceleration zone Bx ~ 1.5 x 10^^ 

A newborn millisecond pulsar spins down due to magnetic dipole radiation in a time 
scale given by = ^-k/^* — ( B'^%n'^ ) ' which for a moment of inertia / = 10^^ g cm^ gives 

~ 4.3 X 10'' s Bi^ ^2 5k- show that the condition that the magnetosphere and 

the disk stresses are in equilibrium at the inner disk edge results a disk mass accretion 
rate Md — ax^ Ic^ (G M^)~^/^QI^^ r~^, where 2 ^ ax > 1 measures the amount of 
magnetic dipole flux that has been pushed by the disk accretion flow to the inner edge of 
the disk (Gosh & Lamb 1978, Shu et al. 1994). Substitution of the previous equations 
yields Md ^ 3 x a2^'^ Bf^Ql^^f, (where a2 — ax/2), which is much larger 

than the Eddington accretion rate Msdd — 7.0 x 1O~-^'''M0 s''^{M^/Mq). However, this 
"super-Eddington" accretion (which is correlated to t^,) will last for a time td, which is 
only a small fraction (fo) of r^. The strong radiation pressure from the accreted material 
will cause most of the infalling material to be ejected from the system, and this will in turn 
cause the accretion rate to rapidly decrease to a value nearly equal to the Eddington limit 
(e.g., Lipunova 1999). Advection dominated inflow-outflow solutions involving supercritical 
accretion onto neutron stars predict a total mass depostion on the star ~ few O.OIM© 
(e.g.. Brown et al. 1999). Thus assuming that a mass M ~ 0.04Mq is accreted during the 
supercritical phase, we find td ~ M/Mj:, ~ 1.3 x 10^ s, and /d ~ 0.03. The acceleration of 
UHECRs in the reconnection zone, on the other hand, will last as long as the supercritical 
accretion. It is well known that the most violent solar flares can live up to several hours 
and the more energetic the longer-lived they are. In our scenario, considering that the 
reconnecting magnetic flelds are several orders of magnitude larger than in the solar corona. 
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it is reasonable to expect that the most violent reconnection events can live at least for 
several days, as required by td. The spectrum evolution of the accelerated UHECRs will 
be, therefore, determined by td = /dt* (see below). 

In order to derive the spectrum of accelerated particles, let us first evaluate 
the rate of magnetic energy that can be extracted from the reconnection region, 
Wb — (Bx/Sn) ^va {4:7rRx Ait!x), where ~ c is the Alfven velocity, and ^ < 1 is a 
factor that determines the amount of magnetic energy released in the reconnection that will 
accelerate the particles. Substituting the previous relations into this equation, one finds 
Wb — 2.6 X 10^^ erg s~^ ^ Bf^ ^zsV According to our model assumptions, in a reconnection 

— * 

burst an electric field arises inductively because of plasma fow across the B lines. The 
corresponding maximum voltage drop is y ~ e Ai?x (vaI c)Bx^Rx — Bx^Rx where e 
is the electric field strength (e.g., Bruhwiller and Zweibel 1992). Once a particle decouples 
from the injected fluid, it will be baUistically accelerated by the electric held and the 
momentum attained will depend on the intensity of the electromagnetic burst and on the 
time and location at which it decouples from the injected fluid motion (e.g., Haswell et al. 
1992). It is out of the scope of this work to quantitatively constrain the burst characteristics 
and thus to deduce the detailed spectrum of particle energies. Clearly, however, the more 
energetic bursts will give rise to a higher average energy for the ejected particles (Haswell 
et al. 1992). In these extreme cases, the induced electric voltage will dominate and the 
particles will be accelerated to an average energy E "^V Ze, which according to the equation 
above and consistently with relation (1) is ~ 10^° eV. The UHECR flux emerging from 
the reconnection site can then be estimated as iV ~ ^ ~ 1.6 x 10^^ ^ B'l^Vt^2 5k ^20 
for particles with energy E ^ 10^° eV, and the particle spectrum N{E) is obtained from 



N = N{E) 



dt — 



N{E) 



I; ^*/fD, or 




) 



1/4 



N{E) ~ 1.7 X 10^^ GeV-^ ^ Z-^^^ B^s'^ E. 



(2) 
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where Eq. (1) has given dQjdE ~ 2.1 x Z'^l^ E~^'^ B:[^'^ {^^fr^) Mh the 



signal made equal in Eq. 1). Eq. (2) above predicts that N(E) (x E ^/'^ = E ^'^j which is a 
flat spectrum in good agreement with observations (e.g., Olinto 2000). 

The particle distribution emerging from the source will not be isotropic. Thanks 
to the magnetic fleld geometry in the reconnection site (see Fig. 1), it will be conflned 
to a ring (above and below the accretion disk) of thickness and height given both by 
~ Ai?x- The distribution will thus be beamed in a solid angle AQ ~ Att{ARx/Rx)'^- 
Let us now, estimate the resulting flux of UHECRS at the Earth. The total number of 
objects formed via AICs in our Galaxy is limited by nucleosynthesis constraints to a very 
small rate ~ (lO"^ — 10~^) yr~^, or in other words, less than 0.1 % of the total galactic 
neutron star population (Fryer et al. 1999). Assuming then that the rate of AICs in the 
Galaxy is taic~^ — 10~^ yr~^, we can evaluate the probability a priori of having UHECRs 
events produced in the Galaxy. The beaming mentioned above will reduce the probability 
of detection of the events of a source by a factor fb ~ (ARx/Rx)"^ — 10~^. Thus the 
probability will be only P ~ fbTAic~^t ~ 2 x 10~^, where t = 20 years accounts for the 
time the UHECR events have been collected in Earth detectors since the operation of the 
first experiments. Since the individual contribution to the observed UHECRs due to AICs 
in our Galaxy is so small we must evaluate the integrated contribution due to AICs from 
all the galaxies located within a volume which is not affected by the GZK effect, i.e., within 
a radius R50 — Rg/50 Mpc. Assuming that each galaxy has essentially the same rate of 
AICs as our Galaxy and taking the standard galaxy distribution no — 0.01 e^^-'^ Mpc~^ 
(Peebles 1993) (with the Hubble parameter defined as Ho = h 100 km s~^ Mpc~^), the 
resulting fiux at £'20 > 1 is F{E) ~ N{E) noTAic'^ Rg, which gives 



( 



ARx/Rx 
0.1 



) 



1/4 



F{E) 



3.3 X lO-^'^CGeV-^cm-^s 




,-3/2 -1 

'20 'TAIC,5 TT-O.Ol -"-50 



(3) 
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where taic,5 — 'Taic ^/lO ^ yr ^, and no.oi = nc/O.Ol Mpc ^. Observed data by the 
AGASA experiment (Takeda et al. 1999) gives a flux at E =10^° eV of F{E) ~ 4 x 10"^° 
Gev~^ cm~^ s~^, so that the efficiency of converting magnetic energy into UHECR should 
be ^ ^ 0.1 in order to reproduce such a signaL 



3. Conclusions and Discussion 

We have discussed the possibihty that the UHECR events observed above the GZK 
hmit are protons accelerated in reconnection sites just above the magnetosphere of very 
young millisecond pulsars originated by accretion induced collapse. AlC-pulsars with 
surface magnetic flelds 10^^ G < ^ 10^^ G and spin periods 1 ms ^ < 60 ms, are 
able to accelerate particles to energies > 10^° eV. These limits can be summarized by the 
condition ;^ 10^^ G (P*/2.5 ms)^^^ (Eq. 1 and Fig. 2). Because the expected rate of AlC 
sources in our Galaxy is very small, the total flux is given by the integrated contribution 
from AIC sources produced in the distribution of galaxies within a volume which is 
unaffected by the GZK cutoff (of radius Rq — 50 Mpc). We flnd that the reconnection 
efficiency factor needs to be ^ ^ 0.1 in order to reproduce the observed flux. This result 
is appeahng because it predicts no correlation of UHECR events with the Galactic plane, 
in agreement with present observations. However, as data collection improves, we should 
expect some sign of correlation with the local distribution of galaxies and the supergalactic 
plane. These predictions can be ultimately tested by coming experiments such as the 
AUGER Observatory which will provide high statistics samples of UHECRs. 

The model predicts a highly super-Eddington accretion mass rate during part of the 
AIC process. In such a regime, one should expect a large mass outflow from the stellar 
surface itself which might alter the coronal conditions near the helmet streamer region. 
However, the strong closed magnetic flelds in the pulsar magnetosphere can inhibit the gas 
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outflow from the steUar surface if the surface temperature does not exceed the value at 
which the radiation energy density {aT"^) equals the magnetic fleld energy density [B^/Sn), 
i.e., ^ 5 X 10^ K bH^; this upper limit is consistent with the predicted values for a 
neutron star that is formed by AIC (Woosley & Baron 1992). 

Let us consider now the energy loss mechanisms that may affect the efficiency of the 
acceleration of the particles to the UHEs. Energy losses by curvature radiation which 
occur when cosmic rays have to stream along magnetic fleld lines with a flnite radius of 
curvature Rc, have a cooling time trc oc R^. In our model, the reconnection takes place 
in the helmet streamer region (see Fig. 1) where i?c — > oo and therefore, trc » ta, 
where ta is the acceleration time for a proton in the reconnection site. [We can use the 
dimensions of the reconnection site derived in §2 to estimate the order of magnitude of 
ta ~ ARx/c ~ 8 X IQ-^s (^^^fr^) ( 245x106cm ) •] The time for synchrotron losses is 
tsyn oc 9"'^, where 9 is the particle pitch angle. Since we require that the protons in the 
reconnection zone are accelerated by an induced electric fleld rather than by a scattering 
process, they are expected to maintain beamlike pitch angles while escaping along the 
magnetic fleld lines, i.e., ^ << 1 and tgyn » ta- 

The accelerated protons may also undergo energy loss by pion and e"^ pair 
production due to interactions with photons from the accretion disk radiation fleld. The 
characteristic distance scales for these processes can be estimated by Ap-y ~ {apjn^)~^ and 
Xpair — {cTpairn^)"^ for photo-piou and pair production, respectively, where is the photon 
number density at the reconnection site, ~ 2.5 x 10~^* cm~^ is the cross section for pion 
production, and apair — 10~^^ cm~^ is the cross section for pair production (e.g., Bednarek 
and Protheroe 1999). In order to estimate n-y, we must determine the luminosity in the 
accretion disk. At super-Eddington accretion rates, the disk is thicker and heat is radially 
advected with matter. The luminosity of such advective supercritical accreting disks is given 
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by Ld ^ [0.6 + 0.7 IniMo/MEdd)] LEdd (e.g., Lipunova 1999), where LEdd - 1-25 x 
erg s"^ (M^/Mq). Using the value obtained in §2 for M^/MEdd — 1-8 x 10^, we find 
Ld — 1.8 X 10^^ erg s~^. Inside the disk, where the optical depth is much larger than 
unity, photons and particles are in thermodynamic equilibrium and the disk radiates like a 
black-body with a temperature :^ 1.1 x lO'^ K {j^^j^jMfr^—r) (ft) {^^j^^) ■ 
At a distance Rx, ^ LD/{nR\e^) (e.g., Sorrell 1987), where ~ ^.^keTo ~ O.SSMeK 
is the mean photon energy, or ~ 8 x lO^^cm-^ (^^^-^^^__^) ' (^^5^) \ 

We note that beyond the disk, the luminosity can be substantially smaller as the matter 
just above the disk surface can be partially opaque to the outgoing radiation. Thus, the 
value above of for the helmet streamer is probably overestimated. Substitution of 
into the equation for Ap-y gives Xp^/Rx ^ 0.22 > 0.1 ( ^^01^^ )' ^ ^^^^ interaction 
distance scale for pion production is larger than the width of the acceleration region. For 
pair production, although the interaction distance Xpair is shorter, the inelasticity of this 
process is so small (inelasticity i^p 2 x 10~^) that the corresponding energy loss rates of 
the UHECR protons are even smaller than in the case of pion production (e.g., Bednarek & 
Protheroe 1999). However, the injected pairs can initiate inverse Compton pair cascades 
in the reconnection region that may saturate an induced electric field if the width of the 
reconnection region is larger than Xgat — 100 (atripTij)'^ , where atrip — 1-3 x 10~^^ cm~^ is 
the cross section for the triple pair production by electron-photon collisions. We find 
that Xsat/Rx ^ 0.4 > 0.1 ( ^"^01^^ )' ^^^^ width of the reconnection site is smaller 
than Xsat- 
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Figure Caption 

Figure 1. Schematic drawing of the magnetic field geometry and the gas accretion fiow in 
the inner disk edge at Rx- UHECRs are accelerated in the magnetic reconnection site at 
the helmet streamer (see text) . The figure also indicates that coronal winds from the star 
and the disk help the magnetocentrifugally driven wind at Rx to open the field lines around 
the helmet streamer (adapted from Shu et al. 1999). 

Figure 2. Allowed zones in the parameter space for UHECR acceleration for E = 10^^ and 
10^^ eV. The vertical fine in each plot indicates the upper hmit on the stellar angular speed 
^-k,max — {G M^l R^y/"^ ~ 1.15 X 10^ s~^. The allowed zone in each plot is above the solid 
line for which ARx/Rx = 1; the dash-dotted line has ARx/Rx = 0.1; and the dashed line 
has ARx/Rx = 0.01. 



